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Abstract 

The nearby (z = 0.03015) cluster of galaxies Abell 2199 was observed by Suzaku 
in X-rays, with five pointings for ~ 20 ks each. From the XIS data, the temperature 
and metal abundance profiles were derived out to ~ 700 kpc (0.4 times virial radius). 
Both these quantities decrease gradually from the center to peripheries by a factor 
of ~ 2, while the oxygen abundance tends to be flat. The temperature within 12' 
(~ 430 kpc) is ~ 4 keV, and the 0.5-10 keV X-ray luminosity integrated up to 30' is 
(2.9 ±0.1) x 10 44 erg s -1 , in agreement with previous XMM-Newton measurements. 
Above this thermal emission, no significant excess was found either in the XIS range 
below ~ 1 keV, or in the HXD-PIN range above ~ 15 keV. The 90%-confidence upper 
limit on the emission measure of an assumed 0.2 keV warm gas is (3.7-7.5) xlO 62 
cm -3 arcmin -2 , which is 3.7-7.6 times tighter than the detection reported with XMM- 
Newton. The 90%-confidence upper limit on the 20-80 keV luminosity of any power 
law component is 1.8 x 10 43 erg s _1 , assuming a photon index of 2.0. Although this 
upper limit does not reject the possible 2.1cr detection by the BeppoSAX PDS, it 
is a factor of 2.1 tighter than that of the PDS if both are considered upper limits. 
The non-detection of the hard excess can be reconciled with the upper limit on diffuse 
radio emission, without invoking the very low magnetic fields (< 0.073/xG) which were 
suggested previously. 

Key words: galaxies: clusters: individual (Abell 2199) — X-rays: galaxies: clus 
ters 
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1. Introduction 



The intra-cluster medium (ICM) in clusters of galaxies, i.e, hot thermal plasmas in col- 
lisional ionization equlibria confined within their gravitational potential, constitutes the most 
dominant form of cosmic baryons that has ever been detected. Although the ICM radiates 
predominantly optically-thin thermal X-ray emission, some clusters have been reported to ex- 
hibit excess signals above the thermal emission at the lowest or highest energy ends of their 
X-ray spectra. One interpretation of the soft excess is emission from Warm Hot Intergalactic 
Medium (WHIM), which has been predicted from cosmological simulations (e.g. Cen & Ostriker 
1999) and is expected to solve the so-called missing baryon problem (e.g Fukugita et al. 1998). 
Another interpretation of the spectral soft and/or hard excess is non-thermal emission from 
accelerated particles in galaxy clusters (e.g. Lieu et al. 1999; Fusco-Femiano et al. 1999): a 
power-low shaped spectrum can exceed the dominant thermal emission at both sufficiently high 
and low energies. This provides a possibility that galaxy clusters are giant accelerators in the 
universe. Thus, the search for soft and hard excess signals in clusters forms a very important 
research subject. However, the instrumental sensitivity in these photon energies has been in- 
sufficient, and hence the existence of these excess components has remained controversial (e.g. 
Bregman & Lloyd-Davies 2006; Nevalainen et al. 2007; Rossetti & Molendi 2004; Fusco-Femiano 
et al. 2007). 

Among those galaxy clusters from which the excess emission has been reported, 
Abell 2199 is particularly interesting, because it has been suspected to exhibit both the soft and 
hard excess components. In fact, Kaastra et al. (1999) reported the detection of both, based 
on the BeppoSAX, Extreme Ultraviolet Explorer, and ROSAT observations. The two excess 
components were simultaneously expressed by a single power law with a photon index ~ 1.8 
which is superposed on the thermal emission spectrum, and interpreted as inverse Compton 
(IC) photons produced when the cosmic microwave background (CMB) photons are scattered 
up by relativistic electrons in the cluster. However, as argued by Kempner & Sarazin (2000), 
the magnetic field in Abell 2199 would have to be unusually weak (< 0.073 /iG) in order for 
the relativistic electrons not to produce synchrotron emission beyond measured upper limits 
of diffuse radio emission. As an alternative interpretation, Kempner & Sarazin (2000) pro- 
posed non-thermal bremsstrahlung by supra-thermal electrons which are being accelerated in 
the cluster. 

In the era of XMM-Newton, the possibility of thermal emission from a warm gas also 
became another popular idea to explain the soft excess in several clusters. Kaastra et al. (2004) 
reported a redshifted O VII emission line from a 0.2 keV warm gas in Abell 2199, along with 
Coma, Abell 1795, Sersic 159-03, MKW 3s, Abell 2052, and Abell 3112, and hence ascribed 
their soft excess to emission from WHIM. At the same time, Chandra has made progress on the 
ICM physics of Abell 2199 especially in the central region, such as an asymmetric temperature 
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distribution in the direction perpendicular to the jets within 30" of the center (Kawano et al. 
2003), and a weak isothermal shock associated with the central active galactic nucleus (Sanders 
& Fabian 2006). 

In order to examine the soft and hard excess phenomena of Abell 2199 with a better 
sensitivity, we observed the object with Suzaku (Mitsuda et al. 2007). The excellent low-energy 
capability of XIS-BI (back illuminated CCD, Koyama et al. 2007) and the low background of the 
silicon PIN detectors in the HXD (Takahashi et al. 2007; Kokubun et al. 2007) have provided 
the best data in the soft and hard X-ray bands, respectively. As a result, we have obtained 
tight upper limits on both excess components, and strengthened the dominance of the thermal 
emission from this cluster. 

In the present paper, quantities quoted from the literature are rescaled to the Hubble 
constant of H = 71 km s -1 Mpc -1 , and errors are given at 90 % confidence range unless other- 
wise stated. 

2. Observation and Data Preparation 

2.1. Observation 

As shown in figure 1 and table 1, we observed Abell 2199 with Suzaku in the beginning of 
October 2006. The observation was conducted in 5 separate pointings ("Center" and "Offsetl" 
through "Offset4"), for ~ 20 ks each, with the XIS field of view (FOV) partially overlapping. 
Throughout the 5 pointings, all the XIS sensors were in normal clock mode without window 
or burst options. These observations became one of the first ones in which spaced-row charge 
injection (SCI) was applied to the XIS detectors (Nakajima et al. 2008). In the present paper, 
we use data of XIS0, XIS1, XIS3 and the HXD-PIN. Although XIS2 was still operational, we 
do not use its data because a database of non X-ray background (NXB) for SCI-ON XIS2 is 
too small, which is available only until the anomaly in XIS 2 on 2009 November 9 1 . 

Between the Offset2 and Offset4 observations, operations of changing the HXD-PIN 
bias voltage were carried out in response to flare-like events in a PIN diode 2 (W10P0). As 
summarized in table 2, data of all the 64 PINs were available throughout the Center and 
Offsetl pointings, and for a part of Offset2 and Offset4. In contrast, only one fourth of the PIN 
diodes were operational in Offset3. Therefore, we do not analyze the Offset3 HXD-PIN data in 
the present paper. For the other observations, we exclude the data of the four PIN diodes in 
W10 unit, and use the remaining 60 channels. Backgrounds and responses excluding the four 
PIN diodes were prepared with a help of the HXD team. 

1 http: / /www. astro. isas.ac.jp/suzaku/doc/suzakumcmo-2007-08.pdf 

2 http:/ /www. astro. isas.jaxa.jp/suzaku/log/hxd/ 
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2.2. XIS data reduction 



The present paper utilizes the XIS and HXD-PIN data, prepared through version 2.0 
pipe-line data processing, and obtained from the Suzaku ftp site 3 . We re-processed un- 
screened XIS event files using Suzaku software version 11 in HEAsoft 6.6, together with 
the calibration datebase (CALDB) released on 2009 January 9. We applied xiscoord, 
xisputpixelquality, xispi, and xistime in this order, and performed the standard 
event screening. Bad pixels were rejected with cleansis, employing the option of chip- 
col=SEGMENT. We selected events with GRADE 0,2,3,4, and 6. Good time intervals 
(GTI) were detemined by criteria as "SAA_HXD==0 && T_SAA_HXD>436 && C0R>6 && ELV>5 
&& DYE_ELV>20 && A0CU_HK_CNT3_NML_P==1 && Sn_DTRATE<3 && ANG_DIST<1 . 5". Details 
of the processing and screening are the same as those described in the Suzaku Data Reduction 
Guide. 4 

The redistribution matrix files (RMFs) of the XIS were produced by xisrmf gen, and 
auxiliary responce files (ARFs) by xissimarfgen (Ishisaki et al. 2007). As an input to the 
ARF generator, we prepared an X-ray surface brightness profile of Abell 2199 using a double-/? 
model (sum of two j3 models; King 1962) of which the parameters were determined through a 
least chi-square fit to the background subtracted and vignetting corrected 0.5-10 keV XMM- 
Newton MOS1 image of the Abell 2199 cluster. The employed parameters are (r c , /3) = (10.3±0.5 
kpc, 0.650 ±0.008) for the narrower component and (51.0 ±0.8 kpc, 0.531 ±0.002) for the wider 
component, where r c is the core radius. Absorption below 2 keV, caused by a carbon-dominated 
contamination material on the XIS optical blocking filters, is included in the ARFs. There, 
differences of the contaminant thickness among the XIS sensors are taken into account, along 
with its radial dependences and secular changes (Koyama et al. 2007). 

Non X-ray background (NXB) of each XIS sensor was created using xisnxbgen, which 
sorts spectra of night earth observations according to the geomagnetic cut-off-rigidity (COR) 
and makes an averaged spectrum weighted by residence times for which Suzaku stayed in 
each COR interval (Tawa et al. 2008). We assumed systematic errors (90% confidence range) 
of 6.0% and 12.5% for the NXB from XIS-FI (front illuminated CCD, Koyama et al. 2007) 
and XIS-BI (Tawa et al. 2008; table 5 (c)), respectively, and added them in quadrature to 
the corrensponding statistical errors. The X-ray background was estimated as described in 
subsection 3.1. 

2.3. HXD-PIN data reduction 

Unscreened event files of HXD-PIN were re-processed by hxdtime, hxdpi, and hxdgrade 
in this order. Using hxdgtigen, we rejected the time intervals when "FIFO Full", "BFSH", or 
"TLMRJCT" happened, and further narrowed the GTI by imposing criteria of "SAA_HXD==0 && 

3 ftp://ftp. darts, isas.jaxa.jp/pub/suzaku/ver2. 0/ 

4 http:/ /ftools. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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T_SAA_HXD>500 && TN_SAA_HXD>180 && C0R>6 && ELV>5 && A0CU_HK_CNT3_NML_P==1 && 
HXD_DTRATE<3 && ANGJDISK1 . 5 && HXD_HV_Wn_CAL>700 && HXD_HV_Tn_CAL>700". 

For the Center, Offset 1, Offset2, and 0ffset4 pointings, we made special HXD-PIN RMFs 
which exclude contributions of the four PINs in W10 unit. This was constructed by summing 
RMFs from the remaining 60 PIN diodes. Since hxdarf gen did not support extended sources 
yet, we created PIN ARFs for each 1' x V pixel by hxdarf gen, assuming the same X-ray surface 
brightness as used for XIS ARFs (subsection 2.2) . The ARFs were then averaged by weighting 
with counts contained in the corresponding pixels. 

As the HXD-PIN NXB, the "tuned" PIN NXB files published by the HXD team 5 were 
used (Fukazawa et al. 2009). In the present paper, we used 2.3% as systematic error (90% 
confidence range) of the PIN NXB. Details of its derivation is described in Appendix 1. The 
PIN CXB was estimated from an observation of a nearby region in the same manner as the 
XIS CXB (subsection 3.2). 

3. Estimation of X-ray Backgrounds 

3.1. XIS X-ray background 

To analyze the XIS data of an extended object in general, we must subtract the Cosmic 
X-ray background (CXB) and Galactic foreground emission (GFE). Since the emission from 
Abell 2199 itself hampers their direct estimation using our own data, we used 73.6 ks Suzaku 
observation of a nearby blank sky region named "High Latitude Diffuse A" (Hereafter HLD; 
Observation ID = 500027010), which is 4°. 03 away from Abell 2199. The difference of their 
Galactic latitudes is only 0°.7, and hence the hydrogen column densities (Dickey & Lockman 
1990, weighted average over 1° cone radius) toward Abell 2199 (0.86 x 10 20 cm" 2 ) and HLD 
(1.02 x 10 20 cm -2 ) are similar. 

After applying the standard processing described in subsection 2.2 to the unscreened 
XIS data of HLD, we removed circular regions (1' in radius) around six point sources that are 
visible in the 0.5-10.0 keV XIS images, and extracted spectra of the remaining region from the 
four XIS detectors. The corresponding detection limit for the point sources is ~ 6 x 10~ 14 erg 
s _1 cm -2 in 2-10 keV. Then, the NXB-subtracted XIS-BI and XIS-FI (averaged over the three 
XIS-FI detectors) spectra in energy ranges of 0.3-5.5 keV and 0.5-5.5 keV, respectively, were 
fitted simultaneously with a CXB+GFE model (see below) using XSPEC12 version 12.4.0ad, 
incorporating the XIS response for a source of uniform brightness. In the analyses described 
below, metal abundances refer to Anders & Grevesse (1989), and the photoelectric absorption 
cross-sections to Balucinska-Church & McCammon (1992) with a new He cross section (Yan et 
al. 1998). 

As the CXB model, we used a power law with a fixed photon index of 1.41 (Kushino et 
5 http://www.astro.isas.jaxa.jp/suzaku/analysis/hxd/pinnxb/tuncd/ 
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al. 2002) and a free normalization. To represent the GFE, we employed model 1 of Henley & 
Shelton (2008), which consists of three apec components (Smith et al. 2001): a non-absorbed 
0.08 keV component, an absorbed 0.11 keV one, and an absorbed 0.27 keV one, representing 
the emission from Local bubble (LB), a cool halo, and a hot halo, respectively. The three 
apec temperatures were all fixed. The hydrogen column density was fixed to the Galatic value 
of 1.02 x 10 20 cm -2 , and the metal abundance and redshift were also fixed to one solar and 
zero, respectively. Relative normalizations of the two halo components were tied as cool halo 
: hot halo = 1 : 0.24 (Henley & Shelton 2008). Thus, the GFE model has only two free 
parameters; normalizations of the LB and halo. As summarized in table 3, this CXB+GFE 
model reproduced the XIS spectra from HLD successfully. The 2.0-10.0 keV surface brightness 
of the CXB component is (6.34±g;ii ± 0.63) x 10~ 8 erg s" 1 cm~ 2 sr" 1 (90% statistical and 
systematic errors), where the systematic error, 10%, refers to Appendix 1 of Nakazawa et 
al. 2009. Using ASCA, Kushino et al. (2002) derived an absolute 2.0-10.0 keV CXB surface 
brightness to be (6.38 ± 0.07 ± 1.05) x 10 -8 erg s _1 cm -2 sr _1 (90% statistical and systematic 
errors). This value becomes (4.2 - 6.0) x 10 -8 erg s _1 cm -2 sr _1 when the surface brightness is 
integrated to our detection limit of 6 x 10~ 14 erg s _1 cm -2 (equation 6 of Kushino et al. 2002). 
Thus, the CXB level obtained from the HLD observation is consistent with the ASCA result. 

3.2. HXD-PIN X-ray background 

We also analyzed the HXD-PIN data of the HLD observation which were processed 
by the same standard procedure as described in subsection 2.3. Ignoring the GFE which is 
negligible at > 3 keV, we fitted the NXB-subtracted 12-40 keV HXD-PIN spectrum with a 
wide-band CXB model. Namely, following the HEAO-1 result (Boldt 1987), we expressed the 
CXB spectral surface brightness as 

(rp \ —0.29 , rp s 

3tev) -H4^v) ergCm " S " Sr " keV "' « 
where So, a normalization, is a free parameter in the fit. When reproducing the CXB signals in 
HXD-PIN, equation 1 was corrected for the 13% systematic difference of cross-normalization 
between the XIS and PIN 6 . The result we obtained, So = (9.7 ±2.0) x 10~ 9 , is consistent with 
the HEAO result of So = 9.0 x 10~ 9 , which can in turn reproduce, within 1%, blank-sky PIN 
spectra which have higher statistics (Fukazawa et al. 2009). The contribution of the six point 
sources to the HXD-PIN spectrum is an order of magnitude lower than the NXB-subtracted 
signal in 12-40 keV. Therefore, the contribution of point sources affects the CXB level only by 
ASo = 0.1 x 10~ 9 , which is much smaller than the statistical error. 

Although we have thus adopted the two separate models to estimate the CXB in the XIS 
and HXD-PIN, the derived two surface brightness results agree within 1% in the connecting 
energy range of 6.0-12.0 keV. Therefore, our two modelings of the CXB are consistent. If we 

6 http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2007-ll.pdf 
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constrain the two CXB models so that the 6.0-12.0 keV surface brightness becomes the same, 
and fit the XIS and HXD-PIN spectra simultaneously, we can constrain the CXB above 10 keV 
better with So = (9.7 ±0.2) x 10~ 9 . This case is shown in figure 2. 

3. 3. Field to field difference of X-ray background 

In the following spectral analysis, we fit the Abell 2199 data simultaneously with the 
HLD data, in order to determine the X-ray background. However, the true X-ray background 
of the Abell 2199 field may be different from that of the HLD field. In order to examine 
the GFE brightness for possible differences between them, we employed X-ray count rates of 
these two fields from the ROSAT All-Sky Survey (RASS) which are available via the NASA's 
HEASARC website 7 . Since the RASS data right on the Abell 2199 field is contaminated, of 
course, with the cluster emission, we obtained a count rate averaged over an annular region 
of which the inner radius is 1°.0 (~ 2.2 Mpc) and the outer radius is 2°.0, centered on the cD 
galaxy NGC 6166. Since the virial radius of Abell 2199 is estimated to be ~ 1.7 Mpc (~ 0°.8) 
from the ICM temperature of 4 keV (Evrard et al. 1996), this annular region can be considered 
to be free from the ICM emission. 

Table 4 summarizes the RASS count rates of the Abell 2199 annular region and the 
HLD field (circular region, 1°.0 in radius). The 3/4 keV (0.47-1.21 keV) and 1.5 keV (0.76- 
2.04 keV) rates are the same between the two fields within statistical errors. Thus, the HLD 
region is confirmed to be a good background estimator for our Abell 2199 observations in these 
energy bands. The 1/4 keV (0.12-0.284 keV) band rates differ by ~ 20%, suggesting a spectral 
difference betweeen the Abell 2199 and HLD fields in the softest energy band. However, the 
LB component, which dominates the spectrum in the 1/4 keV band, is poorly determined with 
the XIS spectra, and the 20% discrepancy is within the statistical error of ~ 40% associated 
with the LB normalization in the following spectral analysis (§ 4.1). 

The CXB surface brightness is known to fluctuate from direction to direction (Ishisaki 
1997). In the case of Suzaku, this fluctuation is estimated to be 10% and 18% for the XIS (full 
FOV) and PIN, respectively (see Appendix of Nakazawa et al. 2009 for details). We take these 
systematic errors into account in the spectral analyses described in section 4. 

4. Data Analysis and Results 

4-1. Thermal emission 

In order to quantify the thermal emission of Abell 2199, we defined seven concentric 
annular regions as shown in figure 1, each with a radial width of 3'. The center is chosen to be 
the X-ray emission centroid at (a,5) = (16 h 28 m 36.9 s , +39°32'53") in J2000.0 coordinates. This 
is 0'.4 offset from the nucleus of the cD galaxy, NGC 6166, but well inside its optical extent of 

http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl 
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~ 2' in diameter. 

In analyzing a given annulus, we extracted XIS-FI (averaged over XIS and XIS 3) and 
XIS-BI spectra from corresponding regions in the five pointings (but excluding those which have 
no intersection), and subtracted the NXB as described in subsection 2.2. For each annulus, 
the utilized observations are given in table 5. Then, the spectra from the same annulus (but 
different fields of view) were read into XSPEC together with the NXB-subtracted HLD spectra 
to constrain the X-ray background. The annular spectra were fitted simultaneously with a 
common vapec model with a free photoelectric absorption, together with the X-ray background 
model described in subsection 3.1, incorporating the responses created for the Abell 2199 ICM 
emission (subsection 2.2) and for a uniform X-ray source, respectively. All the parameters 
of the X-ray background model were tied between the annular and HLD spectra, except for 
the thickness of XIS contamination which depends on the date of observation (Koyama et 
al. 2007). Among the different observations of Abell 2199, all the vapec parameters, except 
the overall normalization, were constrained to be the same. The cross normalization between 
XIS-FI and XIS-BI was set free. Metal abundances, Z meta \, in the vapec model were tied in 
three groups as Zq = Z^ c = ^M g = Zau Zsi = Zs = Zav = Zc&, and Zp c = Zw An example 
of this fitting in the 0'-3' region is shown in figure 3. We also changed the normalization of 
CXB model by ±11% and studied the effect of CXB fluctuation described in subsection 3.3. 
Although the X-ray background model is determined separately for the seven individual regions, 
the normalizations of CXB, LB, and Halo are consistent among the regions within errors (see 
figure 10 in Appendix 3). 

The results of this analysis are summarized in figure 4 (for details, see table 7 in 
Appendix 5). The single-temperature vapec model thus reproduced the spectra successfully, 
and the CXB fluctuation did not affect the results significantly. The temperature and iron 
abundance were determined out to ~ 700 kpc, which is about 0.4 times the virial radius of 
a 4 keV cluster. The temperature gradually decreases toward cluster outskirts from ~ 4 keV 
to ~ 3 keV. The temperature within 12', ~ 4 keV, is consistent with that derived with ASCA 
(4.1 keV; Fukazawa et al. 2004). The iron and silicon abundances also decrease outwards by 
a factor of ~ 2, while the oxygen abundance tends to be flat. The hydrogen column density 
obtained from the fit became (2-3) x 10 20 cm -2 , which is significantly higher than the Galacitic 
value. We examine the cause of this effect in subsection 4.2. The 0.5 - 10 keV X-ray lumi- 
nosity integrated out to 30' is (2.9 ±0.1) x 10 44 erg s _1 , which is consistent with the result of 
XMM-Newton (Snowden et al. 2008). 

Next, we attempted a two-temperature fit using a model expressed as vapec+vapec. 
Metal abundances of the two components were tied together. As shown in figure 5a and 
figure 5f (for details, see table 8 in Appendix 5), the spectra in the innermost region (0'- 
3') were reproduced better with the two-temperature model of 2 keV and 5 keV, while the 
improvement was not significant in the other regions. This is reasonable, because Johnstone et 
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al. (2002) detected, with Chandra, a temperature drop from 5 keV to 2 keV toward the center 
within a central region of radius 100 kpc (2'. 8). The iron and silicon abundance profiles did not 
change significantly between the single-temperature and two-temperature models. However, 
the oxygen abundance in the inntermost region decreased and the profile became flatter. The 
high values of Nr, found with the single-temperature analysis, persisted; this effect cannot be 
regarded as an artifact caused by an inappropriate temperature modeling. 

4-2. Constraints on soft excess emission 

A soft excess may manifest itself as a hydrogen column density Nr that is significantly 
lower than the Galactic HI columnn density toward the Abell 2199 field (weighted- averaged 
value over 1° radius from NGC 6166), 8.60 x 10 19 cm -2 from Dickey & Lockman (1990) or 
8.92 x 10 19 cm -2 from Kalberla et al. (2005). The measured values of Nr in figure 4 and 
figure 5, however, are factor of 2-3 higher than the Galactic value in all regions. While this 
result generally argues against the presence of any flux excess in lower energies, we need to 
understand the nature of the higher absorption before discussing the soft excess issue. There 
are five possibilities to explain this effect: (1) excess absorption within the Abell 2199 system, 
(2) an overestimate of the background which was subtracted, (3) inadequate modeling of the 
ICM emission, (4) presence of high metallicity clouds at the direction of Abell 2199, and 
(5) underestimate of the XIS contaminant thickness which is included in the ARF. The first 
possibility is unlikely, because the measured Nr is high even in the outer regions such as 12'— 15'. 
The second alternative is rejected, because we already incorporated the uncertainties of NXB 
and CXB in the analysis, and because within 6', the background level is more than an order of 
magnitude lower than the signal below 1 keV band. The third is also unrealistic, because the 
values of Nr remained unchanged between the single- and two-temperature modelings. The 
forth is probably unphysical because metal abundances in such clouds would have to be as large 
as ~ 3 solar in order to explain the excess X-ray absorption, while keeping Nr to the Galactic 
HI value. Then, we focus on the last possibility. 

When we obtain Nr separately with XIS-FI and XIS-BI in the three inntermost regions 
where the values of Nr are constrained well (figure 4(b)), the results with XIS-FI become sys- 
tematically higher than those with XIS-BI by ~ 20%. This difference suggests an instrumental 
cause such as (5). Then, we evaluated an excess thickness of contaminant, assuming (5) is 
the cause of the high Nr. We included an additional absorption by the XIS contamination in 
the single temperature model as vapecxphabsxvarabs. Here, the varabs factor represents the 
additional absorption which is added to that already included in the ARF, and can handle sep- 
arately column densities of elements from hydrogen to nickel. The varabs column densities of 
all the elements except carbon and oxygen were fixed to zero, and the relative column-number- 
density ratio of carbon to oxygen was fixed to N c /N = 6 (Koyama et al. 2007). By introducing 
the additional varabs factor, the best-fit model parameters did not change significantly except 
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iVn in phabs. Figure 6 shows the confidence contours between iVn in phabs and the additional 
carbon column density Ac in varabs. As easily expected, these two quantities correlate neg- 
atively with each other, and we cannot determine the true Ah accurately. Fixing Ah at the 
Galactic HI value in figure 6 (dotted magenta lines) yields iVc ~ 5 x 10 17 cm -2 (dashed cyan 
lines in figure 6), which is comparable to an upper limit allowed by its systematic uncertainty 8 . 
Therefore, we may regard the excess absorption observed in figure 4b and figure 5b as due to 
an underestimation of the XIS contaminant by iVc ~5x 10 17 cm -2 , and presume that N-r is 
equal to the Galactic HI value. From figure 6 of the innermost annuli, we can reject any Ah 
value which is significantly lower than the Galactic HI value (dotted magenta lines in figure 6), 
because the extra carbon thickness must be < 5 x 10 17 cm -2 . 

In figure 7, we show data-to-model ratios in the 0.3-1.0 keV range, when the data are 
fitted by the single temperature vapec model with Ah = 8.60 x 10 19 cm -2 and Ac = 5.0 x 10 17 
cm -2 both fixed. The value of iVc is thus set to the allowed maximum, and the combination of 
Nc an d Nn means a slight (~2.5cr) overestimate of the overall absorption specified by the data in 
figure 6 of the innermost annulus. Furthermore, the data are divided by the single-temperature 
model which was determined over a wide band of 0.3-8.0 keV. Therefore, soft excess emission, 
if any, should be seen in these ratio spectra. However, no significant systematic excess is seen 
therein. Then, to obtain an upper limit on the soft excess due to a thermal warm gas, we 
added to the model a 0.2 keV apec component with 0.2 solar metallicity, following the result 
of Kaastra et al. (2004). The derived 90% upper limit on surface brightness of the 0.2 keV 
component is summarized in table 6. In table 6, we also give XIS upper limits when a power 
law with photon index of 2.0 is added to the model. The upper limit on the power law in 
terms of the 0.2-10 keV luminosity becomes 7.1 x 10 43 erg s" 1 (2.5 x 10 43 erg s _1 in 20-80 keV 
range) when integrated to 30'. In subsection 4.3, we independently constrain such a power law 
component using the HXD-PIN data. 

As a cross confirmation, we repeated the same analysis by choosing Ac = and A^h = 
3 x 10 20 cm -2 , the latter suggested by figure 4. However, the ratio spectra remained essentially 
unchanged. 

4-3. Constraints on hard excess emission 

In figure 8, we show a 15-40 keV NXB-subtracted HXD-PIN spectrum summed over 
Center, Offset 1, Offset2, and Offset4. By this summation, the statistical error associated with 
the 15-40 keV HXD-PIN signal reduced to ~ 0.8% of the NXB, which is considerably smaller 
than the systematic NXB error of ±2.3% (90% confidence range) adopted in the present paper. 

After subtracting the NXB, the HXD-PIN signal (black crosses in figure 8) was thus 
detected significantly up to ~ 30 keV. For comparison, we also show in figure 8 spectra of 
the CXB and the thermal emission. Here, the thermal emission was modeled by summing 

8 http://www.astro.isas.jaxa.jp/suzaku/process/caveats/cavcats_xrtxis08.html 
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the single temprature vapec models determined for individual annular regions. As to thermal 
emission outside the XIS FOV which contributes ~ 11%, we assumed that the surface brightness 
obeys the double-/? profile which was described in § 2.2, and that the temperature and metal 
abundances therein is the same as those in the outermost region covered by the XIS. After 
the nominal XIS vs. HXD-PIN cross normalization (Kokubun et al. 2007, see also a Suzaku 
calibration document named suzakumemo-2007-11 9 ), the model normalization of the thermal 
component was multiplied by a factor of 1.13 when calculating its predicted contribution to the 
HXD-PIN data. Using this cross normalization, intensities of the two CXB models in 6.0-12.0 
keV agree within 1% as described in § 3.2. This assures that the systematic errors associated 
with the cross normalization is < 1%, which is much smaller than the statistic erros in the HXD- 
PIN spectrum. Thus, the NXB-subtracted spectrum is consistent, within statistical errors, with 
a sum of the CXB and the extrapolated thermal emission (cyan points in figure 8). Even when 
the NXB is set to the minimum that is allowed by the 90% systematic uncertainty range, the 
HXD-PIN signal does not exceed the (CXB+thermal) emission significantly. Thus, the hard 
X-ray excess was not detected. 

In order to obtain an upper limit on hard excess emission, we lowered the NXB within 
its 90% systematic uncertainty, and also lowered the CXB model within its fluctuation of 18% 
(90%-confidence). The estimation of the CXB fluctuation is described in Appendix 2. Then, we 
fitted the NXB-subtracted spectrum with a power law (with its photon index fixed at 2.0 and 
its normalization left free) plus the fixed CXB and thermal models, and obtained upper limits. 
This very conservative method gave a 90%-confidence upper limit of 1.8 x 10 43 erg s _1 on the 
20-80 keV luminosity of the power-law component. If this component is uniformly distributed 
over a circular region of radius 10', the implied 0.2-10 keV luminosity per solid angle becomes 
1.6 x 10 41 erg s _1 arcmin -2 . This limit is comparable to those obtained with the XIS using the 
softest specctral end. 

5. Discussion 

Using the five pointing observations with Suzaku, we measured the temperature and 
metal abundances of Abell 2199 out to ~ 700 kpc (0.4 times virial radius), and searched the XIS 
and HXD-PIN data for soft and hard excess emissions, respectively, both above the thermal 
component. We however found no significant excess emissions in either energy bands, and 
derived their upper limits. Over a broad energy band of ~ 0.4 keV to ~ 30 keV, the emission 
is dominated by the thermal components with temperatures of a few keV. 

The upper limits on soft excess, expressed as the emission measure per unit solid angle 
from a 0.2 keV warm gas (table 6), are more stringent, typically by more than a factor of 3, 
than the reported positive detection by XMM- Newton (Kaastra et al. 2004), (28 ± 13) x 10 62 

9 http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2007-ll.pdf 
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cm -3 arcmin -2 . One likely cause of this inconsistency between Suzaku and XMM-Newton is 
confusion with solar wind charge exchange (SWCX) emission, which has been observed with 
Chandra (Wargelin et al. 2004), XMM-Newton (Snowden et al. 2004), and Suzaku (Fujimoto 
et al. 2007). In fact, proton flux near the Earth, observed with WIND-SWE 10 , enhanced by a 
factor of 3-4 for a part of the XMM-Newton observatons compared with that of Suzaku. The 
relevant proton light curves are shown in Appendix 4. Similarly, in Suzaku observations of 
outskirt regions of the Coma cluster, Takei et al. (2008) found no significant excess O VII or 
O VIII emissions which were previously reported with XMM-Newton (Finoguenov et al. 2003), 
and concluded that the excess emission is likely due to an enhancement in the SWCX emission 
during the XMM-Newton observations. 

The Suzaku data have constrained the 20-80 keV luminosity of any excess hard emission 
to be < 1.8 x 10 43 erg s _1 , assuming that it has a power-law spectrum of photon index 2.0 and 
employing the most conservative estimates of the CXB. Using the BeppoSAX PDS, Nevalainen 
et al. (2004) detected non-thermal emission from Abell 2199 at 2.1cr level, with 20-80 keV 
luminosity of (1.7 ±1.3) x 10 43 erg s _1 (90%-confidence errors) assuming a power-law spectrum 
of photon index 2.0. This detection by the PDS is not rejected by our result. However, the 
BeppoSAX results, with a significance of only 2.1cr, would be considered as an upper limit, like 
the present result. In this case, our upper limit is a factor of 1.7 more stringent than that of the 
PDS. From a calibration analysis using the Crab Nebula, the 20-80 keV flux derived from the 
PDS is known to be systematically lower by 21% than that derived from HXD-PIN (Nakazawa 
et al. 2009). If we take into account this systematic effect, the difference in the upper limits 
between HXD-PIN and PDS increases to a factor of 2.1. 

Considering the difference between the FOVs of HXD-PIN (34' in FWHM) and PDS 
(1.3° in FWHM), the claimed detection by the PDS becomes consistent with our upper limit 
if the excess hard X-ray emission is much more extended than the HXD-PIN FOV. However, 
the total FOV of HXD-PIN in our pointing observations (~ 43' = 1.5 Mpc radius in FWHM) 
mostly cover the whole cluster region, ~ 0.9 times the Virial radius of Abell 2199 (~ 1.7 Mpc, 
see § 3.3). Therefore, to reconcile the two results, we would have to assume that the emission 
is spatially distributed much beyond the virial radius. 

According to Kempner & Sarazin (2000), a radio synchrotron flux density at 327 MHz, 
S u , expected from Abell 2199, is related to the 10-100 keV IC X-ray flux, S'hxr, as 

Sv = 234 (lO- 11 erg cm" 2 s^ 1 ) (l^c) ( 327 MHz) Jj ' (2) 
assuming the photon index of the hard X-ray emission to be 1.81. Here, B is the cluster 
magnetic field, and v is the observed frequency. The upper limit they quoted on the diffuse 
radio flux, S v < 3.25 Jy at 327 MHz, and the hard X-ray emission detected by BeppoSAX, 
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required a very weak magnetic field (< 0.073 fiG). However, given our non-detection of such 
an excess hard X-ray flux, this is no longer the case. When the HXD-PIN data are used to 
constrain a power-law component with a photon index of 1.81 (instead of 2.0), its 10-100 keV 
luminosity becomes < 1.5 x 10~ n erg cm -2 s -1 . As a result, effectively any strength of magnetic 
field is allowed as shown in figure 9. 

Regardless of the magnetic field strength, the present result constrain the number of 
relativistic electrons in the system, so that their IC emission should not exceed the present 
hard X-ray upper limit. If we assume a power-law spectrum of synchrotron emission with a 
photon index 2.0, the corresponding electron number spectrum has a power-law index 3.0, and 
is written as (Petrosian et al. 2006) 

iV( 7 ) = 2iV tot7min 7" 3 (7>7min), (3) 

where 7 is the Lorentz factor, N tot is the total number of electrons, and 7 m ; n is a lower-cutoff 
in 7 which is conservatively estimated to be ~ 10 3 (Petrosian et al. 2006). These relativistic 
electrons produce IC emission with a 20-80 keV luminosity of 

L HX R = 1.0 x 10 45 (^g) (1 + zf erg s"\ (4) 

where z is the redshift. From equation 4 and our upper limit, we obtain N tot < 1.6 x 10 63 . Then, 
the integrated electron kinetic energy K e for 7 > 7 m ; n in Abell 2199 becomes 

K e = N tot m c c 2 (2 lmin - 1) < 2.6 x 10 60 (^jf) erg, (5) 

where m e is the electron mass. Since electrons with 10 3 ^ 7 ^ 10 4 , which contribute to the 20- 
80 keV IC emission, have a liftime of ~ 10 8-9 yr in typical ICM conditions (Petrosian 2001), 
the acceleration luminosity in Abell 2199 can be constrained as < 8.4 x 10 44 (7mm/10 3 ) erg 
s _1 . This upper limit is rather loose, and is three times higher than the 0.5-10 keV luminosity 
of the thermal emission. If we assume that the magnetic field of Abell 2199 is ~ 1 fiG as 
suggested by Faraday rotation measurements (Fe & Owen 1994), the radio upper limit gives 
the dominant constraint on the relativistic electrons, with the predicted Shxr two orders of 
magnitude below the present upper limit (figure 9). In this case, the constraint of acceleration 
luminosity becomes < 10 43 (7 m in/10 3 ) erg s^ 1 which is about an order of magnitude lower than 
the 0.5-10 keV luminosity. 

The total mass of Abell 2199 is 1.3 x 10 14 M Q when integrated to 0.8 Mpc (Fukazawa et 
al. 2004). If Abell 2199 was formed via a major merger between two clusters with similar masses 
M (~ 7 x 10 13 M Q each) at a relative speed of v ~ 3000 km s _1 , the kinetic energy deposited 
onto the merged system becomes 

terser ~ ~Mv 2 = 6 x 10 63 ( ^ | ( r v) erg. (6) 

mcrger 2 \7x 1O 13 M ; V 3000 km s^ 1 / & y 1 

Since shocks produced in cluster mergers are thought to last typically r ~ 10 9 yr (Takizawa 
2000), the energy input rate to non-thermal electrons becomes E- m ~ 1 x 10 46 (//0.05) erg s -1 , 
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where / is the fraction of energies given to them. Values of E- m = io 46-47 erg s _1 are also 
expected from theoretical models of merging clusters such as Coma (e.g. Takizawa 2004). 
Then, our upper limit on the acceleration luminosity is an order of magnitude lower than E- m , 
assuming / ~ 0.05. We hence infer that Abell 2199 have been free from major merger events for 
more than ~ 10 9 yr, which is a typical time scale on which the IC emission decays (Takizawa 
2004). 

MK acknowledges support from a grant based on the Special Postdoctoral Researchers 
Program of RIKEN. The present work is supported in part by the Grant-in-Aid for Scientific 
Research (S), No. 18104004. 

Appendix 1. Estimation of NXB systematic error 

la - systematic erros of PIN 15-40 keV band, using earth occultation data and E0 102-72 
observations which has been regularly scheduled for the XIS calibration, is 2.3% (Fukazawa et 
al. 2009). Since the statistical error is 1.8% (ler), the systematic uncertainty of the PIN NXB 
model in 15-40 keV band is estimated to ±2.3% in 90% confidence region. 

Appendix 2. Estimation of CXB systematic error 

The CXB has uncertainty because the number of unresolved faint sources in FOV statis- 
tically fluctuate. This CXB fluctuation scales as Q~ ' 5 S®' 25 , where Q e is the FOV and S c is the 
detection threshold flux for point sources (Ishisaki 1997). By scaling from the CXB fluctuation 
of HEAO-I, 2.8% for Q c = 15.8 deg 2 and S c = 8x 10~ n erg s" 1 cm" 2 (Shafer 1983), that of the 
HXD is estimated to be 18%, assuming Q c = 0.32 deg 2 and S c = 8x 10~ 12 erg s" 1 cm -2 in the 
10-40 keV band. This corresponds to 1.3% of the NXB level. 

Appendix 3. Normalization of the background model 

The background model (GFE+CXB) are determined separately for the seven annular 
regions. We show the normalizations of background components in figure 10. The normaliza- 
tions of background models are consistent among the regions within errors (90% confidence 
range) . 

Appendix 4. Proton flux of Abell 2199 observations 

figure 11 is proton flux of WIND-SWE 11 , when XMM-Newton and Suzaku observed 
Abell 2199. These plots are created by multiplying proton speed and proton density of which 
data are public in the WIND-SWE web site. 

11 http : / / web . mit.edu/space/ www /wind . html 
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Appendix 5. Fitting results of the seven annular regions 
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Fig. 1. A 0.5-10 keV mosaic XIS image from the Center, Offsctl, Offset2, Offset3, and Offset4 obser- 
vations, smoothed with a gaussian of a — 10" and corrected for vignetting. Boundaries of annular regions 
used for the spectral analysis are shown in green circles. See the electronic version of the paper for a colour 
figure. 
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Fig. 2. NXB-subtractcd XIS-FI (black; 0.5-5.5 keV), XIS-BI (red; 0.3-5.5 keV), and HXD-PIN (green; 
12-30 keV) spectra of the HLD observation, fitted with the X-ray background model described in the text. 
The CXB is modeled by a power-law in the XIS band, and equation 1 in the HXD-PIN region, with their 
surface brightness constrained to match over 6.0-12.0 keV. A model for the GFE is included in the XIS 
range. 
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Fig. 3. NXB-subtracted XIS spectra of Abell 2199 (0' - 3' region) and the HLD. Spectra from Center 
XIS-FI (black; 0.4-8.0 keV), Center XIS-BI (red; 0.3-8.0 keV), Offset 1 X1S-F1 (green; 0.4-8.0 keV), and 
Offsctl X1S-B1 (blue; 0.3-8.0 keV) are fitted with a common vapec model and the X-ray background 
model. The XIS-FI (cyan; 0.5-5.5 keV) and XIS-BI (magenta; 0.3-5.5 keV) spectra from the HLD field 
are fitted simultaneously with the same X-ray background model. Further details are described in the 
text. 
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Fig. 4. (a) The temperature profile of Abell 2199 determined by the single temperature vapec model 
(black). Those when the CXB level is set maximum and minimun are also shown in red and green, 
respectively (b) The same as panel (a), but for the hydrogen column density. The Galactic value is 
shown as a horizontal dotted line, (c)-(e) Abundances of Fe, Si, and O, respectively, with the same color 
specifications as panel (a), (f) Reduced chi-squared of the fit. The degree of freedom is also shown in the 
parentheses. See the electronic version of the paper for colour figures. 
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(a) Temperature profile (b) Hydrogen column density profile (c) Fe abundance profile 
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Fig. 5. The same as figure 4, but for the two temperature vapec fits. Blue data points in panel (a) 
indicate the lower-temperature component. In the 3'-6', 9'-12', and 12'-15' regions, they were fixed at 
0.5 times the corresponding best-fit higher temperature because they were not constrained. Magenta and 
cyan data points are the lower-temperature component when the CXB level is set maximam and minimum, 
respectively. See the electronic version of the paper for colour figures. 
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(a) 0'-3' (b) 3'-6' (c) 6'-9' (d) 9'-12' 
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Fig. 6. Confidence contours calculated for the seven regions between the hydrogen column density Nn, 
and the extra carbon contaminant thickness Nq which is added to that already included in the nominal 
XIS ARFs. Confidence levels are ler (black), 2a (red), 3er (green), and 4a (blue). The Galactic HI colunmn 
density of the Abell 2199 field, 8.60 x 10 19 cm -2 (Dickey & Lockman 1990), is shown in a dotted magenta 
line. A typical systematic uncertainty of the contaminant thickness in the ARF, 5 x 10 17 cm~ 2 , is also 
shown in a dashed cyan line. See the electronic version of the paper for colour figures. 
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Fig. 7. Data-to-model spectral ratios in the 0.3-1.0 keV range, when the data in each annulus are fit- 
ted with a single temperature vapec model with fixed hydrogen (iVu = 8.92 x 10 19 cm -2 ) and carbon 
(N c = 5.0 x 10 17 cm' 2 ) column densities. From top to bottom, the ratio spectra of the 0'~3', 3'— 6', 6'-9', 
9'-12', 12'— 15', 15'-18', and 18'-21' regions are shown. Colors specify different observations. See the 
electronic version of the paper for a colour figure. 
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Fig. 8. NXB-subtracted 15-40 keV PIN spectra (black) summed over the Center, Offset 1, Offsct2, and 
Offsct4 observations. The same spectrum when the NXB is lowered within the 90% confidence limit is also 
plotted in red. The CXB model (green), the extrapolated thermal emission (blue), and their sum (cyan) 
are also plotted for comparison. 
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Fig. 9. Relations between the 10-100 keV hard X-ray flux and the 327 MHz flux density of Abell 2199, 
produced when the same population of relativistic electrons with index 2.6 interact with the cosmic 
microwave background photons and magnetic fields of various strengths, respectively. The upper limits 
on the radio (Kempner & Sarazin 2000) and hard X-ray (this work) signals are shown in solid lines, and 
the gray region satisfies these limits. 
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(c) HALO normalization 
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Fig. 10. (a) normalization of the CXB determined by the single-temperature vapec model (black). Those 
when the CXB level is set maximum and minimum are also shown in red and green, respectively, (b) The 
same as panel (a), but for normalization of the LB. (c) The same as panel (a), but for normalization of the 
Halo, (d) The same as panel (a), but for the the two-temperature vapec model, (e) The same as panel (b), 
but for the the two-temperature vapec model, (f) The same as panel (c), but for the the two-temperature 
vapec model. See the electronic version of the paper for colour figures. 
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Fig. 11. Proton flux (speed times density) measured by the WIND-SWE when XMM-Newton and Suzaku 
observed Abell 2199. 
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Table 1. Suzaku Observation log of Abell 2199. 



Position (Obs. ID) 


Start* 


End* 


XIS Exposure 1 ' 


PIN exposure^ 


Center (801056010) 


2006/Oct/01 07:08:00 


2006/Oct/01 21:31:18 


18.5 


19.3 


Offsetl (801057010) 


2006/Oct/01 21:32:15 


2006/Oct/02 10:43:19 


23.5 


23.7 


Offset2 (801058010) 


2006/Oct/03 12:41:00 


2006/Oct/04 03:47:19 


19.0 


10.8 


Offset3 (801059010) 


2006/Oct/04 03:48:52 


2006/Oct/04 16:17:19 


18.9 




Offset4 (801060010) 


2006/Oct/04 16:19:16 


2006/Oct/05 07:18:24 


23.5 


15.3 



*: Time is shown in UT. 



': Effective exposure in units of ks, obtained after the event selection. 



Table 2. HXD-PIN HV operations during the observation*. 



Position 


HVP0 


HVP1 


HVP2 


HVP3 


Center 


400 


500 


500 


500 


Offsetl 


400 


500 


500 


500 


Offset2 


400 


500 -► 


500 -> 


500 -» 


Offset3 


400 











Offset4 


400 


-> 400 


-> 500 


-> 500 



*: The numbers indicate high voltages in units of volts. 



Table 3. Results of model fittings to XIS spectra of the HLD observation. 





t(cxb) 


&(CXB)t 


Sx(GFE)* 


X 2 /dof (prob.§) 


1.02 x 10 20 (fix) 


1.412 (fix) 


6-35±g;i| x lO" 8 


i.96iS;i? x io- 8 


315.1/285 (0.11) 



*: Hydrogen column density in units of cm" . 

The 2.0-10.0 keV CXB surface brightness in units of erg s _1 cm -2 sr 
*: The 0.3-1.0 keV GFE surface brightness in units of erg s -1 cm -2 sr~ 
=*: Null hypothesis probability. 
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Table 4. RASS averaged count rates around Abell 2199 and the HLD field. 



Band* 


1.0°-2.0 ot 


HLD* 


1/4 keV 


1454.7 ±8.2 


1745.1 ±17.4 


3/4 keV 


144.0 ±2.8 


143.1 ±5.9 


1.5 keV 


133.4 ±2.7 


139.4±4.6 



*: RASS 1/4 keV, 3/4 keV, and 1.5 keV bands correspond to 0.12-0.284 keV, 0.47-1.21 keV, and 0.76-2.04 keV, 
respectively. 

t; RASS count rate in an annulus around Abell 2199, in units of 10 -6 counts s -1 arcmin -2 . 

RASS count rate averaged over the HLD field (1° radius) in units of 10 -6 counts s -1 arcmin -2 . 



Table 5. Observations used in the XIS annular spectrum fitting. 



Re 


gion 


Observation 


0'- 


3' 


Center, Offsctl 


3' 


6' 


Center, Offsctl, Offset2 


6' 


9' 


Center, Offsctl, Offset2 


9' 


12' 


Center, Offsetl, Offset2, Offset4 


12' 


-15' 


Offsetl, Offset2, Offset3, Offset4 


15' 


-18' 


Offsetl, Offset2, Offsct3, Offset4 


18' 


* 


Offset3, Offset4 



Table 6. Upper limit on the soft excess. 



Region 


0'-3' 


3'-6' 


6'-9' 


9'-12' 


12'-15' 


15'-18' 


18'-* 


Warm gas* 


7.5 


7.3 


4.0 


6.5 


3.7 


6.0 


5.3 


Power law* 


3.8 


1.2 


0.7 


0.5 


0.3 


0.4 


0.3 



*: Surface brightness of 0.2 keV gas, expressed as emission measure per unit solid angle in units of 10 62 cm 
arcmin 

0.2-10 keV luminiosity per unit solid angle of a power law with photon index of 2.0, in units of 10 41 erg s 
arcmin -2 . 
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Table 7. 



Fitting result of the single-temperature vapec model for the seven annular regions. 



Rcj 


jion 


kT* 








7 % 




7 % 

z o 


X 2 /D.O.F. 


CXB level is default 


0' 


-3' 


o ni+0.04 
J ' al -0.04 


2 


71+0.21 

' ' x -0.20 


u - oz -0.02 


66+ 007 





70+0.10 

'°-0.09 


2687.5/2189 


3' 


6' 


4 2fi+ 05 

^■ ZD -0.05 


2 


^n+0.24 
•^°-0.23 


n c.q+0.02 


n 07+0.08 

u -°' -0.08 





c-r+0.12 
oo -0.12 


2014.8/1823 


6' 


9' 


4 in+ - 06 


2 


or+0.50 

z,J -0.26 


U - OD -0.03 


o.34i8;ll 





rc+0.18 
oo -0.12 


1295.1/1067 


9' 


-12' 


o nn+ 012 
o.»u_ 10 


1 


oo+0. 46 
■ 88 -0.57 


n oq+0.05 
u -°°-0.04 


0.3818:1? 





65+ 27 

oo -0.23 


790.5/732 


12' 


-15' 


o ^7+0.13 
J - ' -0.21 


2 


nfi + 1.03 
UO -0.56 


0.3818:85 


20+ - 27 

U.ZU_ 20 





cq+0.44 
UO -0.26 


625.1/552 


15' 


-18' 


o R1+0.15 
°- ul -0.38 





sn+1-41 
ou -0.50 


o.3o«:o9 


< 0.24 





79+0.56 
' ,z -0.47 


539.3/456 


18' 


-21' 


o 1Q+0.36 
°- iu -0.34 





on+0.95 
SU -0.80 




< 0.33 





r ? +0.87 

°'-0.53 


448.5/406 


CXB level is minimum 


0' 


-3' 


o ni+0.04 
°- yi -0.04 


2 


71+O.2I 
' ' x -0.20 


u - oz -0.02 


o.66i8:8 7 7 





70+O.IO 
'°-0.09 


2688.0/2189 


3' 


C 


4 26+ - 05 

^■ zo -0.05 


2 


4 o+0.24 
^ J -0.23 


n QQ+0-02 
u - oa -0 .02 


n 07+0.08 

u -°' -0.08 





rr+0.12 

OO -0.12 


2015.6/1823 


6' 


9' 


4 10+ - 07 


2 


cm +0.39 
°^-0.34 


u -°°-0.03 


33 +0 11 

U - OJ -0.11 





cc+0.16 
oo -0.13 


1295.6/1067 


9' 


-12' 


o 04+O.H 


1 


oo+0. 50 
' 8z -0.51 


n oq+0.05 
u -°°-0.05 


39+ - 18 

U - OJ -0.17 





fi c+0.26 
OO -0.22 


790.6/732 


12' 


-15' 


3.56ig;S 


2 


qo+0.64 
' oo -0.88 


n 07+0. 08 
u -° ' -0.06 


1Q+ - 27 

u - iy -0.19 





cc+0.40 
DO -0.19 


624.8/552 


15' 


-18' 


3.68i^? 





on+1-53 
JU -0.30 


u - ou -0.09 


< 0.24 





fiq +0.60 
DJ -0.47 


539.3/456 


18' 


-21' 


07+0.38 

°- z '-0.33 





71 +0.95 
' ± -0.71 


U - JU -0.14 


< 0.33 





°^-0.58 


446.1/406 


CXB level is maximum 


0' 


-3' 


ni+0.04 
°- yl -0.04 


2 


71+0.21 
' ' x -0.20 


u - oz -0.02 


66+ 007 

u - DD -0.07 





70+O.IO 
'°-0.09 


2686.9/2189 


3' 


6' 


4 26+ - 05 

^■ ZD -0.05 


2 


4 o+0.24 


QO+0.02 
u - oy -0.02 


n 07+O.O8 
u -°' -0.08 





c-r+0.12 
oo -0.12 


2013.9/1823 


6' 


9' 


4.0718$ 


2 


qo+0.36 
oo -0.36 


36+ 03 

u - JD -0.03 


33+° 11 





cc+0.16 
oo -0.14 


1293.9/1067 


9' 


-12' 


3-87±g:g 


1 


oq+0.52 
' oo -0.49 


r, qq + 0.05 

u -°°-0.05 


40+ 017 

u -* u -0.09 





65+ 27 

oo -0.21 


790.3/732 


12' 


-15' 


40+0.13 

J -^°-0.19 


2 


1 1+1.03 

■^-o.ss 


n 07+0. 08 

u -°'-0.08 


20+ - 28 

u - zu -0.20 





64 +0.42 
'°<±-0.27 


625.6/552 


15' 


-18' 


4 r+0.19 





4C + 1.59 

'^°-0.45 


o.28i8:g 


< 0.26 





79+0.47 
' z -0.50 


539.4/456 


18' 


-21' 


2 89+ - 37 

Z - Oa -0.38 





oo + l. 04 
■ S8 -0.88 


o.25±8:^ 


< 0.32 





r-q+0.80 
do -0.53 


452.2/406 


*: Temperature 


in units of keV. 














Hydrogen column density 


in " 


units of 10 2C 


1 —? 

cm . 











Metal abundance in solar unit of Anders & Grevesse (1989). 
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Table 8. 



Fitting result of the two-temperature vapec model for the seven annular regions. 



Region 


kT* 


n H Z Fo 


7 t 
Z Si 


7 X 

z o 


X 2 /D.O.F. 


CXB level is default 



0' 


3' 


4 


•791^/2 


04 +o.n 

u ^-0.19 


2 


74+0.21 

' ,4: -0.21 


0M»% 


66+ 06 





40+0.09 

'^°-0.09 


2550. 


.8/2187 


3' 


6' 


5. 


orl^A 


65 (fix) 


2 


4fi +0.24 
/±o -0.23 




n 07+0. 08 

U - J ' -0.08 





45+0.11 
■ 40 -o.io 


1990 


,2/1823 


6' 


9' 


4 


■6oi° : 9 3 l/i 


7,4+0.94 
■ '+-0.39 


2 


49+0.38 
'^ z -0.37 


n oq+0.03 
U - JO -0.03 


o.34±° :i? 





OO+0.17 

JO -0.16 


1278 


.6/1065 


9' 


12' 


4. 




.45 (fix) 


1 


oq+o.51 

oa -0.50 


q + 005 


0-40l: ;i« 





cq+0.24 
do -0.23 


788. 


.1/732 


12' 


-15' 


4. 


05±g;l|/2. 


18 (fix) 


2 


q K +0.76 
JD -0.75 


0.331°;°, 


91+0-25 
u - zl -0.21 





40+0-36 

/±J -0.32 


624 


,0/552 


15' 


-18' 


3 




oq+0.43 


1 


09+0. 98 
■ Oz -0.97 


0.30t° ;il 


<0.48 





r 7 +0.59 

■°'-0.54 


529 


,1/454 


18' 


-21' 


3 


■34l°;«/o 


C4+O.I6 
°*-0.21 


1 


■56-1.23 


°- 37 tai7 


<0.51 





r 7+0. 94 
°'-0.56 


439. 


,6/403 



CXB level is minimum 



0' 


3' 


4 


.79^/2.04+"^ 


2 


7R+0.21 
'°-0.21 


0.481°;^ 


fifi+°° 6 





4 9+0.09 
• 4S -0.09 


2572 


.8/2187 


3' 


6' 


5. 


07l°;°*/2.64 (fix) 


2 


4O+0.24 
' 4S -0.24 


0.38lS:g 


n 07+0. 08 
u - J '-0.08 





45+o.n 


2019 


.8/1823 


6' 


9' 


4 


.631^/1.761°^ 


2 


40+0. 38 
' 4S -0.37 


n oq+0.03 
U - JO -0.03 


S3+ 012 

u -°°-0.11 





qO + 0.17 

JO -0.16 


1314 


.9/1065 


9' 


-12' 


4. 


651^3/2.41 (fix) 


2 


00+0-51 
uu -0.50 


S0+ 05 

u - ou -0.05 


39+ 017 

u - oy -0.17 





c.q+0.24 
do -0.23 


794 


.5/732 


12' 


-15' 


4. 


13l°;g/2.16 (fix) 


2 


r 7 +0.76 
■°'-0.74 


0.33l°;° 7 


1Q+ - 25 
u - ly -0.19 





4O+0.35 
/±a -0.32 


634 


.0/552 


15' 


-18' 


3 


.791^/0.901°- 


1 


R9+0.97 
Oz -0.97 


0.30l° i 4 4 


<0.44 





55+0.48 
oo -0.53 


534 


.4/454 


18' 


-21' 


3 


^1+0-42/0 54+O.IO 
■ oi -0.39/ u -°^-0.20 


1 


79 + 1.39 
' ' z -0.55 


QQ+0-23 
U - OM -0.18 


<0.49 





fifi +0.87 


455 


.7/403 



CXB level is maximum 



0' 


-3' 


4.791^/2.0410^ 


2 


79+0.21 

' ' z -0.21 


o^s±8:8S 


fifi+°° 6 


n 40+0.09 


2571 


,3/2187 


3' 


-6' 


5.061^9/2.66 (fix) 


2 


,0+0.24 
*°-0.23 


0.38l° ;° 2 


o.38lS:g 


0.45l° ;i 


2000 


,3/1823 


6' 


-9' 


4 55+0.94/-, 71 +o.97 

^• oo -0.33/ L - ' 1 -0.37 


2 


qc+0.38 
Od -0.37 


n oq+0.03 
U - JO -0.03 


0.341°;- 


0.381°;^ 


1277 


,4/1065 


9' 


-12' 


4.55l°;i/2.49 (fix) 


1 


75+0.51 
'°-0.50 


qi+0.05 
u - J1 -0.05 


41 + - 18 

u - 41 -0.17 


5T K) - 24 

U.OO_ 23 


820 


,6/732 


12' 


-15' 


3.90l°;g/2.30 (fix) 


2 


10 +0.77 
1U -0.75 


0.33l° ;° 7 


n 9Q+0.25 
U ' ZO -0.23 


0.511^ 


646 


,2/552 


15' 


-18' 


3.6ll°;i/0.89l°;« 





qr+0.99 
ad -0.95 


0.30l° ;il 


<0.49 


0.58l° J° 


551, 


,0/454 


18' 


-21' 


3.20lgS/0.54lg : g 


1 


99 + 1.33 

■ zz -1.22 


0.331°;- 


<0.53 


0.54l°; 9 5 I 


443, 


,1/403 



*: Temperatures of (higher component) / (lower component) in units of kcV. 

T; Hydrogen column density in units of 10 20 cm" 2 . 

*: Metal abundance in solar unit of Anders & Grevesse (1989). 



